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Abstract: The spectra of masses and decay constants for non-strange meson resonances 
in the energy range 0-2.5 GeV is analyzed. It is known from meson phenomenology that 
for given quantum numbers these spectra approximately follow linear trajectories with a 
universal slope. These facts can be understood in terms of an effective string description for 
QCD. For light meson states the trajectories deviate noticeably from the linear behavior. 
We investigate the possible corrections to the linear trajectories by matching two-point 
correlators of quark currents to the Operator Product Expansion (OPE). We find that the 
allowed modifications to the linear Regge behavior must decrease rapidly with the principal 
quantum number. After fitting the lightest states in each channel and certain low-energy 
constants the whole spectrum for meson masses and residues is obtained in a satisfactory 
agreement with phenomenology. We briefly speculate on possible implications for the QCD 
effective string. 
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1. Introduction 

Hadron phenomenology tells us that the masses squared of mesons with given quantum 
numbers lie approximately on equally-spaced linear trajectories^ (see, for example, the 
modern reviews [1,2]). This is a strong indication that QCD admits an effective string 
description, as this type of spectrum is characteristic e.g. of the bosonic^ string [3]. Long 
ago it was noticed [5] that if one accepts this behavior, by including an infinite number 
of resonances it is possible to reproduce the parton-model logarithm present in two-point 
correlators. 

^From now on, we shall term this behavior "linear trajectory ansatz" (LTA). 

^This string model is of course not consistent, a more realistic amplitude is the Lovelace-Shapiro ampli- 
tude but unfortunately it cannot be derived from any known string and as we shall see is also incompatible 
with QCD anyway. See e.g. [4] for a discussion on this point. 
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This observation means that the exchange of an infinite number of vector mesons 
may be dual to the perturbative QCD continuum made of quarks and gluons. This quark- 
hadron duality was explicitly verified in QCD in 1+1 dimensions in the large- A'^c limit [6,7]. 
At present, the QCD community understand by quark-hadron duality several, somewhat 
different, things. The most extended assertion is that perturbation theory can be used to 
calculate a certain smeared hadronic cross section (so-called global duality) [8]. We use 
the term in a rather lax sense and understand by quark-hadron duality any matching of 
hadronic and partonic degrees of freedom. 

In the context of large- A^'c QCD, quark-hadron duality for mesons has received a lot 
of attention [9]- [23]. In [9] the parameters of linear vector mass spectrum were derived 
within the framework of finite-energy sum rules. The axial-vector channel was also con- 
sidered there. Later, this approach was extended to the pseudoscalar [10] and scalar [11] 
cases. In [13] the problem of matching this infinite set of meson states to the Operator 
Product Expansion (OPE) [14] was considered, with my(n) being a general function. In 
the review [15] possible deviations from quark-hadron duality were investigated (including 
meson resonance widths, i.e. away from the large-AT^ limit.) In [16] the linear trajectory 
ansatz in the vector and axial channels were used to calculate certain physical observ- 
ables. Subsequently, in [17]- [23] various aspects of quark-hadron duality, chiral symmetry 
restoration (CSR), and matching to the OPE have been discussed. 

Since, after all, it is well known that the linear trajectory ansatz fails to describe the 
lightest states in each channel (particularly so in the pseudoscalar channel), some attempts 
to go beyond the linear trajectory ansatz were studied in [18, 19], in the latter paper using 
a semi-classical string analysis. The ansatze proposed there, however, seem to lead to a 
discrepancy with the OPE that will be explained later. 

In the present work we want to re-take these issues. Our interest has, to some extent, 
been triggered by the recent controversy among the authors of references [17] and [21], but 
we shall also report here on some older results of us [22]. We shall consider radial Regge 
trajectories for meson resonances with the same quantum numbers (so-called "radial exci- 
tations") in the vector (V), axial- vector {A), scalar (S), and pseudoscalar (P) correlation 
functions. If one believes in the existence of a QCD effective string — and we do — the 
slope of trajectories must be the same in all channels because this slope is proportional to 
the string tension that depends on gluodynamics only (string universality). In this respect 
we tend to agree with the approach of [17], but we find other aspects of this work to be 
more questionable. 

We propose a systematic method to take into account possible corrections to the linear 
trajectories in the V, A, S, P channels. Our ansatz for the meson mass spectra, once sub- 
stituted into two-point correlators, is matched to the OPE, which gives a set of constraints 
on meson mass parameters. As we shall see, the OPE dictates a very particular class of 
corrections to the linear trajectory ansatz. These corrections are not anecdotic: they turn 
out to be essential to get good phenomenology and by taking them into account one does 
not need to introduce an artificial separation between low-lying and excited radial excita- 
tions in the different channels like in [16,17,20,21]. Finally, the results are fitted to the 
experimental data [1,2,24]. 
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The paper is organized as follows. In section 2 we remind the general idea of how 
to construct the asymptotic sum rules for the case of infinite number of resonances. In 
particular, the status and the necessity of non-linear corrections to the linear trajectory 
ansatz is thoroughly examined in sections 3 and 4 for the case of V, A mesons. In section 5 
this analysis is extended to the S, P channels. Section 6 contains the details of our fits and 
discussion of the obtained results. In section 7 we outline our conclusions. 



2. Current-current correlators and scheme of sum rules 

Let us consider the two-point correlators of V, A, S, P quark currents in the large- A^c limit 
of QCD. On the one hand, by virtue of confinement they arc saturated by an infinite set 
of narrow meson resonances, that is, they can be represented by the sum of narrow meson 
resonances with given quantum numbers J and masses mj{n), n = 0, 1, . . ., 



U'iQ^) = [ /xexp(iQx)(gr9(x)grg(0))p,a„ax = ^ n2^^^''l ^ + + D(Q^ (2.1) 

J = S,P,V,A; r = i, 75, 7;i, 7/^75; ^o, Di = const. (2.2) 

The last two terms both in the 5, P and in the V,A channels represent a perturbative 

contribution, with Dq and Di being contact terms required to eliminate divergences in 
the ultraviolet. On the other hand, their high-energy asymptotics is provided [14, 25] by 
perturbation theory and the OPE. In the chiral limit^ 

n (<3 ) = 7-0 1 + — 1^7^ + TTTz — T^A ir^"^^vr ' (2-3) 



47r2 V ' TT y Q2 ■ i27r 9 ' Q6 

,P^r^2^_ 3 / , llas\^2, , {{G^f) , 14 {qqf 



where we have defined 

n):^^(Q2) = (_<5^^q2 ^ g^Q^) n^'^(Q2) . (2.7) 

The vacuum saturation for the 4-fermion condensate [14] in the \aige-Nc limit has been 
used. In Eqs. (2.3)-(2.6) ((G^^)^) and the {qq) represent the gluon and the quark conden- 
sate, respectively. The normalization scale ji arises after additive renormalization of the 



^In [26] asymptotic terms governed by the dimension- two gluon condensate ("gluon mass") were 
introduced. This dimension-two condensate cannot be produced by a local gauge-invariant operator [27,28]. 
On the other hand, this condensate brings httle influence on the fits of meson parameters (see the discussion 
in Appendix A). We put here this condensate to zero. 



-3- 



infinite sums in Eq. (2.1). as is the QCD coupling constant taken at the Chiral Symmetry 
Breaking (CSB) scale ~ 1 GeV. 

As it can be seen from Eqs. (2.3)-(2.6) differences of correlators of opposite parity 
decrease very rapidly at large momenta 

n^Q') - U^iQ') = -STra.M! + O (^) , (2.8) 

This signifies the chiral symmetry restoration (CSR) at high energies, with the quark 
condensate (qq) being an order parameter of the chiral symmetry breaking in QCD. 

Expanding the right-hand side of Eq. (2.1) in inverse powers of and comparing 
with Eqs. (2.3)-(2.6) one obtains several asymptotic sum rules for each channel (we do not 
account for condensates of dimension eight and higher). If one deals with a finite number 
of resonances the procedure is straightforward (for the CSR sum rules [29,30]) . In the case 
of infinite number of resonances and given an ansatz m?{n) one uses the Euler-Maclaurin 
summation formula or, alternatively, the ^/j-function (both methods are equivalent). It 
should be noticed here that, generally speaking, applying a summation procedure for a 
divergent sum and expanding in are not commutative operations. A regulator is in 

general called for. 

The sum rules, obtained in this way, deliver certain constraints on the parameters of 
meson mass spectra. This information can be extracted and compared with phenomenology, 
which is the subject of the next sections. 

In this work we shall ignore the anomalous dimensions of the operators and the run- 
ning of the coupling constant, and therefore we do not take into account the logarithmic 
corrections that they may introduce in the sum rules. We shall come back to this issue in 
a separate paper. 

3. Linear trajectories for vector and axial-vector resonances 

In the V,A case the residues Zv,A{n) in Eq. (2.1) have the structure 

Zy,A(n) = 2F^^^(n), (3.1) 

where the quantities Fv,Ain) represent the decay constants parameterizing the matrix 
elements of the corresponding currents 

(Olje^m(O)l^(e, /c,n)) = .J-^eFv{n)mv{n)e^', (3.2) 

(27r)2 

{0\Ai^{0)\A{e,k,n)) = -^-^eFA(n)mAin)e''. (3.3) 
(27r)2 

where k is the meson momentum, is a polarization vector, and e is electron charge. The 
quantity Fy (n) can be obtained from the decay of vector mesons 
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where a is the fine structure constant. The constants Fji{n) are related to the widths of 
some radiative decays. Unfortunately, this sector is poorly known, except for the reactions 
including the ground axial-vector state (ai-meson). We have found two examples [31,32] 

GlmlFl ( ml\f 2^2 \ ^ ^ 



and 



2 \ 3 

(3.6) 



where Gp is the Fermi constant, is the mass of the r-lepton, and is the pion decay 
constant. 

Let us now introduce the linear trajectory ansatz for the meson mass spectrum 

= Mv,A + av,An ; FyA{n) = const = Fy j^; n = 0, 1, 2, . . . , (3.7) 

with constant intercepts My^ and slopes av,A- In order to do the sums over all resonances 
we shall use the Euler-Maclaurin summation formula in Eq. (2.1) 

^ /-TV 1 D 

^/(n)= / /(x)dx + -[/(0) + /(Ar)] + -f[/'(7V)-/'(0)]- 

n=0 

§[/'" (iV) - f"'m + ... + (-1)' (2fT^[/^''^'^(^) - /^''^'HO)] + . . . , (3.8) 

where Bi = 1/6 , B2 = 1/30 ,.. . are Bernoulli numbers. We shall expand the final expres- 
sion in powers of Q~'^ and after comparing with the OPE's (2.3) and (2.4) one gets the set 
of asymptotic sum rules. 

At this point we have to introduce a regulator^ since the series is manifestly not 
absolutely convergent (the general term behaves like ~ 1/n). We shall cut-off the infinite 
sum by including only the N first terms in the sum. This naturally cuts off the restricted 
value of momenta to ^cut-v A = ^VA + o,v,aNv,a — clv,aNv,a- As we identify parity- 
doublers between the V- and A-channels and do not intent the cutoff to break chiral 
symmetry, the cutoff numbers Ny and Na must be taken equal, Ny = Na- On the other 
hand, the appropriate cutoffs in the OPE also coincide in the V, ^-channels in a chirally 
symmetric regularization^ , so Acut;V = Acut;A (see also similar arguments in [18]). Therefore 
for a natural implementation of chiral symmetry the slopes of radial Regge trajectories must 
be the same, ay = Qa- To compare with OPE (2.3) and (2.4) one has to make an additive 
renormalization by means of subtraction of infinite constants D^'^. In order not to break 



* Another way to improve the convergence is to differentiate with respect to Q^. This recipe will be used 
below for the matching to QCD. 

^Note that a momentum cut-off is questionable in the OPE expansions on account of gauge invariance , 
but perfectly gauge invariant in the sum over resonances (2.1). So one is actually forced to use different reg- 
ulators on both sides, but after the renormalization with the help of constants D^'"^ the leading logarithmic 
term in is in any case unambiguous. 
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the chiral symmetry the subtraction must be equal for parity doublers, 
conclusion is then 

4 Nc _ 2F^ _ 2Fl 
3 167r^ av a-A 

Wc remark that the QCD string model implies ay = ga = a as well. (This relation 
av = aA is not fulfilled in [16, 20] where such analysis for the LTA was performed.) Then 
one concludes that Fv = Fa- Note that the limit iV ^ oo, or equivalently, — > oo is 
assumed. This allows us to neglect terms of the form Q^/A^u^ or i^'^/^cut- 

However, as we have already mentioned, it is not possible to extract further conse- 
quences from the individual V and A asymptotic sum rules. Sub-leading terms are some- 
what ambiguous. Fortunately, the annoying logarithm is absent in differences of correlation 
functions. Namely 

n''(gvn-(g=)^2|(^^-^^) (3.10) 

This difference (and the equivalent for S, P) show a very fast CSR. The leading contribution 
is of 0{l/Q^): all perturbative and purely gluonic contributions cancel in the difference. 
For individual sum rules we shall perform the appropriate derivatives in before matching 
OPE and radial Regge trajectories. 

The above fast convergence in correlator differences was used among others by [17], 
where the following generalization of the Weinberg sum rules [33] was proposed 

00 

J2 K(n)m2j(n) - Flin)mfin)) = C«, i = 0,l, (3.11) 

n=0 

^(0) ^ ^(1) ^ 0. 

Wc emphasize that one sums over chiral pairs in the difference (3.11). This becomes very 
important at the moment that one needs to cut-off the sums with a finite value of A^. If 
the cut-off is placed in such a way that the "chiral partner" of a given resonance, that is 
included, is left out, chiral symmetry will be explicitly broken by the regulator. It is thus 
vital to properly identify "chiral partners". This is not an issue in the V,A cases, but it 
will be so as we shall see in the S, P channels. 

From the second sum rule (i = 1) it immediately follows, given that Fa = Fy and 
a A = o,v, that Ma = My. This automatically implies if the believe in the LTA that the 
spectrum of QCD is degenerate in the vector and axial-vector channels, something that it is 
obviously not true (even after correcting for the quark mass differences, not included in the 
present analysis). Furthermore, the first sum rule (i = 0), implies = 0, something again 
manifestly incompatible with experimental evidence. Since fj^ can be interpreted as the 
order parameter of chiral symmetry, what we are seeing here is that the LTA corresponds to 
a string model that has many common features with QCD, but an essential one is missing. 



= D^. The 
(3.9) 
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namely it would correspond to a string model without chiral symmetry breaking . Note 
that convergence of the infinite sums in Eq. (3.11) also requires to impose the equality of 
slopes ay = a^, intercepts My = Ma, and decay constants Fy = Fa- 

Of course one may escape the contradiction between the fact that the LTA leads to 
/tt = and the real world by declaring that the ground states (perhaps even the first 
few resonances) in each channel, e.g. the p and ai mesons, are isolated resonances. This 
pattern roughly meets the existing spectroscopy of V, A mesons, but it does not reproduce 
the physical values of condensates [16]. However, such an ansatz satisfies the CSR (3.11) 
at high energies. Obviously this way out is quite ad-hoc and cannot be justified on any 
QCD-inspired string model. 

In the present paper we want to argue that the ansatze proposed in [16, 17] can be 
improved by introducing certain systematic corrections to the linear trajectories (3.7). The 
problem here is that arbitrary ansatze for m^{n) and F^(n) result in appearance of terms 
which arc absent in the standard OPE, namely, terms with a fractional power of and 

terms of the kind Q'"^^ ^^7^1^ • (Recall that we do not consider anomalous dimensions.) 



4. Sum rules and deviation from the LTA in the V, A channels 

Our goal is to construct a class of radial Regge trajectories that does not lead to the 
unwanted terms. Namely, our ansatz must reproduce the parton-model logarithm and 
contain only inverse powers of large momentum squared besides the logarithm. Recall 
that, using (3.8), 

Clearly the logarithm can be only produced by the integral. In order to generate the 
logarithm we require that 



/ 



^'W"^ +D„ = ClJ9l±^]^o(']. (4.2) 



Q^ + m^{x) " \ iJ,^ J \Q 



here C is a constant from OPE (3) and (4), Dq is a subtraction constant, and /x is a 
normalization scale. 

Differentiating Eq. (4.2) with respect to x, one can see that the two requirements are 
satisfied only if 

F\x) = tix)^:^^^^^^, t{x) = C + At(x) , (4.3) 
dx 



®The Lovelace- Shapiro amplitude, after correcting the intercept has the proper Adler zero, required 
by chiral symmetry. However, it still has regularly spaced levels, something that as we see is manifestly 
incompatible with QCD asymptotics. It cannot be the correct QCD string. The interested reader may want 
to see [4] for an attempt to formulate consistent string propagation in a background where chiral symmetry 
has been broken. 
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where At (x) is a decreasing function to be defined. If we do not consider anomalous 
dimensions the function At (x) does not have to induce any logarithms and other non- 
polynomial in terms in the integral (4.2). Thus, the direct expansion of the integral 



/ 

Jo 



"»'(^) At (x) d (m2(x)) _ ^ / 1 ^ 



where 

tk= At {x) {-m\x)f d {m\x)) , (4.4) 

Jo 

must exist at any order. 

Convergence of the integrals in (4.4) for any k is possible if the function At {x) falls 
off as an exponential, or faster. Therefore, since one does not have yet any dynamical 
arguments let us take the simplest possibility, namely 

At{x) = Ape-^'"'' , Bf>0, (4.5) 

with constant Ap and Bp. Given an ansatz for m^{n), Eq. (4.3) together with (4.5) 
provides the condition of consistency with the OPE. It is clear that the LTA ansatz (3.7) 
is a simple particular case of Eq. (4.3). Any improvement of the ansatz (3.7) has to meet 
Eq. (4.3) as well. We remark that the corrections to linear trajectories, proposed in [18,19], 
do not satisfy this requirement. Thus, these ansatzc cannot be matched to the OPE. 
The sum rules (3.11) can actually be generalized to values of i greater than 1 

oo 

J2 {FUn)m''^{n) - Fl(n)mi^(n)) = CW, i = 0, 1, . . . . (4.6) 

n=0 

Now the C*^*-* contain a contribution from higher dimensional condensates. For the absolute 
convergence of Eq. (4.6) at a given i one has to have 

1 



F^{n) - Fl{n) ^ a > i, (4.7) 



m^(ra) - m\{n) - ^, /3 > i. (4.8) 
Let us now discuss corrections to Eq. (3.7) 

m^^^(n) = M^_^ + an + Sv,A{n) . (4.9) 

Note that the convergence for any i requires My = Ma = M and that the contribution 
from the function 6v,A{'n) to the difference (4.8) must decrease at least exponentially. 

Collecting our observations, we propose the simplest ansatz satisfying all previous 
requirements, keeping of course in mind that this is just one of the possibilities 

mlr^^in) = M^ + an + A^'^e'^"''', (4.10) 
FlAn) = a(c + A^^\-^^A, (4.11) 



87r2 TT 



C = ^(l + ^), Bm>0, Bp>0, (4.12) 
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with certain constants A^l"^, Bm,F to be fitted. We do not know the underlying dynamics 
responsible for the appearance of those exponential corrections. But it seems for us reason- 
able to suppose that, for masses, this dynamics is governed mostly by gluons and thereby 
does not depend on flavor. Thus, we keep the exponent 5^ the same for all channels. For 
the same reason we regard Bp as independent of parity. We note also that in Eq. (4.3) it 
is enough to retain only the linear in n part of m?{n). 

In what follows we adopt a perturbative approach and retain in the sum rules only 
terms linear in the exponentially small corrections. Products of exponentials exceed the pre- 
cision of our simple (one-exponential) parameterization: these products are regarded as of 
the order of next-to-next-to-leading corrections to Eqs. (4.10) and (4.11). One should keep 
in mind the precision of our approach and thus, when comparing the OPE and the sums 
of resonances, we only match the perturbative logarithms and the first non-perturbative 
contribution 0{1/Q'^) in the V,A channels separately, and the leading (non-perturbative) 
correction 0{1/Q^) in the difference (Q"^) — Il'^{Q'^) which represents a true order pa- 
rameter of the chiral symmetry breaking in the chiral limit. 

In order to avoid of irrelevant infinite constant let us consider the first derivative of 
V, ^-correlators. Introducing the notation for the linear part, 

fh^{n) = M'^ + an, (4.13) 

and making use of (3.1) and (4.10)-(4.12) we have (the indices V, A, are dropped for brevity) 

dm') _ of- a{C + Ape-^--) 
'^Q^ ^(Q2 + m2(n) + A^e-^-'^)'' 

For the ground states the exponential corrections are not small (in general, about 50%). 
We will not apply our perturbative approach for these mesons. Separating out the ground 
states and retaining in the remainder only the parts linear in exponential corrections, one 
has for (4.14) 

ldn(Q2) aiC + Ap) 



2 dQ2 ^Q2j^M^ + Am) 



\2 



^ll(g2 + m2(n))2 (g2 + ^2(„))2 (Q2 + ^2(„))3j \ " ) 

In the r.h.s. of Eq. (4.15) one has three sums. The first one represents the first derivative 
of V'-function, which has the standard asymptotic expansion at large . In the second and 
third sums one may permute the summation over n and expansion in due to absolute 
convergence in any order of this expansion. The final expressions for the sum rules are 
presented in Appendices B and C. (see Eqs. (B.3)-(B.7)). 

The solution of these equations will be considered after discussion of the scalar sec- 
tor. In addition, the so-called D-wave vector mesons should be taken into account, see 
Appendix D. 

As a further evidence for the need of nonlinear corrections let us consider the chiral 
symmetry restoration limit (gg) and /tt 0. As follows from Eqs. (B.3),(B.4), in 
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the case of LTA the corresponding solution = 5*^ gives a negative value for the gluon 
condensate in Eqs. (B.5), (B.6). This implies instability of vacuum energy [34] and, hence, 
the appearance of tachyons. In that sense the LTA is inconsistent with rather general 
principles, even in the chiral symmetry restoration limit. A consistent ansatz can be only 
non-linear. 



5. Sum rules for scaleo" and pseudoscalao* resonances 

The scalar case differs from the vector one and the derivation presented in the previous 
section needs to be revisited. Due to the transversality (2.7), the V, A correlators required 
only one subtraction in Eq. (2.1). The S, P correlators require two subtractions. Accepting 
the definition^ 

Zs,p{n) = 2G%p{n)m'sp{n), 
the analog of (4.2) in the 5, P case takes the form 



(5.1) 



G'^(x)m^(x) dx 

q2 _|_ jYi'^(^x) 



+ Dq + DiQ^ 



-CQ^ In 



7 



+ 



(5.2) 



where C > 0, and Dq,Di are subtraction constants. Let us make the following rearrange- 
ment. The sum over resonances can be rewritten as 



n(Q2) = 2^ 



G^(n)m^(n) 

Q2 _|_ j7^2^j^^ 



+ Do + DiQ^ 



2G^(n) 

Q2 _|_ jYi2(^^'j 



(5.3) 



The first sum in the right-hand side of (5.3) represents an infinite constant and it must be 
renormalized in a chirally symmetric way 



^2G2(n)+Do = ^o, 



(5.4) 



where Dq is different for the S and P channels. Then we have 



Q'G\x)dx l„o2_^o2lr, 



Q'^ + m^{x) 
7 



+ ... . 



(5.5) 



Repeating the discussion in the previous section, which has led to Eq. (4.3), one obtains 
in the scalar case 

(5.6) 



Glp{n) =a(c + ^g'^e-^^'^) , Be > 0, 



'^This formula does not work for the 7r-meson if we include the pion in the radial Regge trajectory. Thus 
for the corresponding residue we just accept its value in the current algebra: Z.,^ = 2^^. As we anyway 
separate out the lowest state, this does not affect the subsequent analysis. Alternatively one could extend 
Eq. (5.1) with the help of a constant shift for all resonances in Zs,p{n) by a quantity Z-,r- However this 
modification induces an unacceptably large value of dimension-two condensate in the OPE (see discussion 
in Appendix A). In accordance to known theoretical and phenomenological estimations [26,27,36] we put 
this condensate to zero. 
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where 

Following the same arguments concerning the parameter B^, as in the previous section, 
we propose the following parameterization of S, P masses 

m|^p(n) = M^ + an + A^^e"-^"^". (5.8) 

Now let us derive the asymptotic sum rules. In order to avoid irrelevant infinite terms 
we consider the second derivative of S, P-correlators. Making use of (4.13), (5.1), (5.6), 
and (5.8) we can write 

{dQ^? ^0 (Q2 + m2(n) + Ame-^rnnf ' 

where the notation rfi^^n) = _|_ an for the linear part of the S,P trajectory has been 
introduced. Separating out the first term in the sum and retaining only the leading order 
in exponential corrections, one obtains 

ld2n(Q2) aiC + Ac) 



4 (dg2)2 (g2+M2 + A„)^ 

aCrh^in) aw?{n)AGe-^Gn 2>aCm^{n)A,„^ 



n=l 

In what follows the procedure is the same as for vector channels. 

At this point one has to decide which particles are chiral partners in order to guarantee 
chirally symmetric results. One may think of two possibilities: a) the vr-meson belongs to 
the radial Rcggc trajectory being the parity-odd partner of the lightest scalar meson, 
b) it does not belong to the radial Regge trajectory being an isolated Goldstone chiral 
particle and therefore the lightest scalar meson is the parity-even partner of the 7r'(1300) 
meson. In case a) the effective low-energy theory would be the linear cr-model, and in case 
b) it corresponds to the nonlinear a-model. We are going to check both variants. The 
corresponding sum rules are presented in Appendix B. 



6. Details of fits and results 



In order to compare the sum of the resonances with the OPE we have to choose appropriate 
values for the inputs of the latter. They are the condensates {qq) and ((G^^/)^), and as- 
On the resonance side we take as input the pion decay constant 7,^, the pion pole residue 
Z-K = 2 fl (from current algebra), and the slope a. The numerical values for all these 
parameters are taken at about the CSB scale ~ 1 GeV and presented in Appendix F. In 
particular, the value^ of /tt ~ 103 MeV is certainly different from its low-energy limit [35] 

®This number can be associated to the trade mark of "Brandy de Jerez 103, Osborne" 
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and corresponds to the matching of the OPE and the sums of resonances at the latter scale 
(see [37] and references therein). 

Recall that in order to keep the discussion in simple terms and also to stay as close as 
possible to the chiral limit we consider here only non-strange mesons. So far we have also 
omitted any reference to the isospin degrees of freedom, but we have now to specify the 
particular pattern in which chiral symmetry is restored in each channel. 

In the V, A channels we study isovector states, where the restoration of U symme- 
try should be manifest in the large Nc limit. We shall therefore analyze the spectra of p 
and ai mesons and their radially excited relatives. In this case five asymptotic sum rules 
contain seven spectral parameters to be fitted: M, Bjn,F, Am^, and A^'^. Thus, we need 
two additional physical inputs, which are chosen to be the masses of the ground states p 
and ai. Actually, we followed the best fit principle when taking the ai mass of about 1200 
MeV in concordance with the masses of first excited vector and axial-vector states. We 
notice also that the ai mass is very sensitive to the variation of /jr, say, its decrease in 3 
MeV (up to 100 MeV) leads to diminishing the ai mass to 1180 MeV. 

Let us comment on some numbers of Table 2 in Appendix F. A possible candidate for 
the state my (2) could be the /9(1900) meson [24]. However, the width of this state is one 
order of magnitude smaller than the widths of other vector states and it could correspond 
to a hybrid state (hybrids and glueballs do not lie on the large- iVc radial Regge trajectory 
for quarkonia). Thus, we do not include this meson on the radial Regge trajectory. 

Unfortunately, we cannot compare the residues of excited states with experiment. The 
possible exception is the vector channel due to Eq. (3.4). For example, our ansatz predicts: 
rp(i45o)^e+e- = 2.9 KeV. The relevant widths Vy_^f.+^- are poorly known and not listed 
in the Particle Data [24]. However, one may compare the results with other independent 
model estimations. The corresponding numbers we have found in [31]: ^ p{i45Q)-ie+ e- — O-^ 
KeV, and in [38]: rp(;^45o)^e+e- = 3.5 KcV, 

Let us discuss now the D-wave vector mesons. Introducing these states entails the 
appearance of three new parameters in the asymptotic sum rules: Md, Ad, and Bd. 
The first one can be fixed by the mass of p(1700)-meson. One can also fix, say -Fy(O) 
and Fd{0). The average of existing estimates of electromagnetic width for /9(1700)-meson 
(^,0(1700)^6+6- = 0-1 KcV [31] and ^p(i7oo)^e+e- = 2.7 KeV [38]) presuppose^ a rather 
large value for Fd(0), of the order 60 MeV. However, we could not find any reasonable 
solution for such large values of residue of ground D-wave vector meson. Our estimations 
showed that -Fd(O) (i.e.. Ad) must be smaller by about two orders of magnitude than 
the value of Fd{0) obtained from the previous analysis. Then the contribution of D-wave 
vector mesons to the physical quantities turns out to be below the accuracy of large- AT^ 
counting. 

^Strictly speaking, Eq. (3.4) can be applied only to the S-wave vector mesons because of locality of vector 
current in (3.2). However, in the relativistic theory the situation changes [31] and transitions (3.2) for the 
D-wave vector states and (3.3) for the ax;ial-vector mesons (which are P-wave states) become possible. 
Moreover, in the large- A/c limit Eqs. (3.2) and (3.3) are well defined even in the relativistic theory since the 
resonances are narrow. Wc use (3.4) for the D-wave vector mesons as a rough estimate of coupling these 
states to the e"''e~-annihilation. 
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In the S, P channel the experimental situation is a lot more confusing. In the large 
Nc limit one might equally look for the restoration of U{1)a symmetry relating isotriplet 
pseudoscalar and scalar states, i.e. pions (/ = 1) and oq (/ = 1) mesons and their radial 
excitations. On the other hand, one could examine the restoration SU{2)l x SU{2)ji 
symmetry relating pions (/ = 1) and /o (/ = 0) in ground and excited states. Both 
scenarios arc compatible if one assumes that there is a degeneracy between oq and /q 
states. To some extent the latter isoscalar [39] can be related to /o(980) from the PDG 
data [24]. But for the former isotriplet there is no firm identification with ao(980) from 
the PDG data [24]. Rather the unitarized fits of pion scattering data [39,40] seem to 
indicate the dynamical origin of ao(980) as a meson bound state. Thus in order to deal 
with well established states in the S, P channels in the present paper wc have restricted 
our analysis to asymtotic sum rules relating the isovector pseudoscalar channel (i.e. pions) 
to the isoscalar meson channel (i.e. /o states). One should of course bear in mind that 
including strange quarks and moving from SU (2) to SU (3) could bring in some relevant 
changes in the scalar sector. 

In this sector we have considered two possibilities, namely the cases a) and b) at the 
end of the previous section. Case a) is labeled "vr-in". By this we mean that the pion 
as well as the lightest scalar meson arc taken to be Regge states. We have here seven 
spectral parameters: M, i?m.,G; i and , which arc subject to three asymptotic 
sum rules. As was pointed above, Bm is universal for all channels. The three additional 
conditions are chosen to be: ms{G) = 1 GeV , mp{0) = 0, and mp(l) = 1.3 GeV. This 
presupposes that the pion and the above scalar are chiral partners. Popular linear cr-models 
often require a much lighter scalar with a mass around 600 MeV. However our fit favors 
a heavier scalar quarkonium to provide a realistic value of Lg, in agreement with the 
conclusions of [39,40]. The numerical results are presented in Appendix F. 

In case b) the 7r-meson is assumed not to belong to the radial Regge trajectory being 
a Goldstone boson that is actually decoupled from the radial Regge trajectory in the CSR 
limit. This case is labeled "7r-out". According to our fits (the lightest scalar state is again 
assumed to have a mass of about 1 GeV), the main qualitative difference between the 
linear and nonlinear cases is that the first one predicts a scalar meson with a mass of about 
1.44 GeV (an iso-singlet resonance with a mass of 1.5 GeV exists but it is widely believed 
to be a glueball), while the second one does not. 

We emphasize again that our analysis in the S,P channels has been performed for 
SU (2) multiplets including iso-triplet pesudoscalar and iso-singlet scalar mesons. We have 
also neglected current quark masses adopting the chiral limit. The situation may be dras- 
tically different for isotriplet scalars and for SU{3) multiplets. In particular, the ao(980) 
(isotriplet) mesons may well be dynamical resonances which decouple in the large Nc limit 
whereas the mesons ao(1450) may be dominantly^^ ground quarkonium states [42]. Their 
rather large masses could arise owing to a strong mixing between a lighter dynamical scalar 
ao and a heavier quarkonium state with masses of order 1.2 GeV [41,42]. We plan to 

^°When lowering the input mass of the scalar meson to 600 MeV one gets Ls = 2.8 • 10~^, much higher 
than the phenomenological estimate Ls = (0.8 ± 0.3) ■ 10"''' [35] 

'^'^We are grateful to the referee who has drawn our attention to this possibility. 
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investigate this scenario with the help of asymptotic sum rules in a forthcoming paper, 
taking also into account light quark mass effects. 

The quantities Gp{n) are related to the corresponding weak decay constants Fp{n) 
through the relation 

Fp(n) = ^Jlh^. 

mp[n) 

In particular, accepting the average value niq for the current masses of u- and d-quarks 
to be equal to 6 MeV, we have the following estimates for the -F,r(i300): Fp{l) = 1.7 MeV 
("-TT-in") and Fp(0) = 1.9 MeV ("7r-out"). Both values are consistent with some previous 
theoretical predictions [43,44], smaller than the estimates from the Finite- Energy Sum 
Rules [10] and larger than the predictions from the non-local quark model [45]. 

It should be also noticed that, allowing for a 10%-accuracy, the first two chiral pairs of 
resonances saturate the chiral constant Lio introduced in [46] almost completely But 
this is not the case for the quantity Atott^ here, for our fits, one needs to retain about seven 
pairs. The reason is that the value of this quantity is very sensitive to the violation of the 
asymptotic sum rules with gluon condensate (strictly speaking, the difference of these sum 
rules) : if one likes to calculate the contribution of the first N pairs of V, A resonances to 
Am^, one has then to deal with exactly this number of resonances in the asymptotic sum 
rules from the very beginning [47]. 

7. Summary 

In the present work wc have considered the matching of the vector, axial-vector, scalar, 
and pseudoscalar meson mass spectra m^(n) {n is the radial quantum number) from Regge 
theory with universal slope to the Operator Product Expansion of quark currents. The 
analysis has been carried out for the light non-strange mesons in the large-iVc and chiral 
limits. Let us summarize the important lessons that we gained from our analysis. 

• The matching to the OPE cannot be achieved by a simple linear parameterization of 
the mass spectrum, the linear trajectory ansatz. 

• The convergence of the generalized Weinberg sum rules requires the universality of 
slopes and intercepts for parity conjugated trajectories. 

• There must exist deviations from the linear trajectory ansatz triggered by chiral 
symmetry breaking. These deviations must decrease at least exponentially with n. 

• There are also deviations from constant residues (decay constants) F'^{n) (or for the 
quantities G^(n) in the scalar case). The analytic structure of OPE imposes again 
an exponential decrease on these deviations (or faster). 

• For heavy states, the D-wave vector mesons have to decouple from asymptotic STim 
rules. This fact implies the exponential (or faster) decreasing the corresponding decay 
constants F'^{n). 

'^^The analytical formulas for the chiral constants Lg, L\o as well as for Am,r are given in Appendix E 
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• Our results seem to exclude a light o"(600) particle as a quarkonium state and rather 
favor the non-linear realization of chiral symmetry with the lightest scalar of mass 
~ 1 GeV, its chiral partner being the 7r'(1300). 

• As a consequence of our approach the quantities Lq,Liq and Am^ are obtained, in 
satisfactory agreement with the phenomenology. 

Unfortunately, the underlying dynamics, which generates the non-linear contributions 
to the spectra of meson masses and residues, is not well known. We can only say that these 
deviations from the string picture seem to parameterize the chiral symmetry breaking in 
QCD [4] and, hence, must be proportional to powers of the chiral condensate {qq) [23]. 
Developing a theory of these non-linear contributions is an interesting task for future. 

Another interesting problem to be examined concerns the SU{3) extension with inclu- 
sion of current quark masses into consideration and a possible resolution of the oq meson 
ambiguity [42]. Finally, one could get some insight by overlapping the present, direct- 
resonance approach with spectral density methods [48]. 
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A. Note on the role of dimension-two condensate 

Let us examine the possibility of having a dimension-two gluon condensate and show its 
unimportance for fitting meson parameters, at least, in the laigc-Nc limit. The dimension- 
two gluon condensate < ("tachyonic gluon mass") was introduced in [26]. This 
dimension-two condensate cannot be characterized by a local gauge-invariant operator: 
for speculations concerning origin, measurement, and physical meaning of this condensate 
see [26-28]. 

In [26, 27] the following relevant modification of the OPE for quark currents was pro- 
posed 

An^(Q2) = An^(Q2) = • ^ (A.i) 

An^(Q2) = An^(g2) = -^A^ In^. (A.2) 

There is no problem to derive such terms from the resonance sums in the V,A channels: 
their introduction is compatible with the ansatze (4.3), (4.10) (sec the corresponding sum 
rules in the next Appendix, Eqs. (B.3), (B.4)). However, in order to reproduce the pertinent 
asymptotical terms in the scalar channels one has to essentially modify the linear part of 
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the ansatz (5.1) for the residues of the S, P resonances. We consider two types of possible 
contributions: an asymptoticaUy constant shift and one affecting only the physical pion 
residue (n = 0), 

Zs,p{n) Zlp{n) = 2G\p{n)mlp{n) + + Sn,o. (A.3) 

For this ansatz the appropriate resonance sums (5.3) generate the contribution 

An^.-c.GoMn(^^i^)+..., (A.4) 

which saturates the dimension-two asymptotics (A. 2) if 

G'o = -||A^ (A.5) 

Then some changes are to be done in the sum rules of the next Appendix. 

On the other hand the modification (A.3) evidently affects the pion pole residue Zt^ 
(at n = 0) 



Z^ = 2^ = Z^ + aG^. (A.6) 



If we put Z,r = and the pion is put on the radial Regge trajectory the required value 
of the "gluon mass" should be ~ — 2GeV^. However this value is, at least, one order 
of magnitude higher that any known theoretical estimations [27], = —(0.2 -i- 0.5)GeV^, 
and phenomenological bounds (from the analysis of r-lepton decay) [36], A^ = —(0.05 ± 
0.08)GeV2. 

We conclude that for realistic values of the dimension-two condensate A^ its contribu- 
tion to the residues Eqs. (A.3), (A.6) is negligible. Furthermore one can check that it is 
certainly less than 5% for meson masses and decay constants. Still a good open question 
is about what is a physical observable which is sufficiently sensitive to its presence. 

B. Sum rules 

After all summations in (4.15) one arrives at the expansion 

k=l 

Substituting the expressions for and comparing (B.l) with OPE (2.3) and (2.4) we 
obtain the following asymptotic sum rules in the V, ^-channels. 

At i/g2 

'l + — ) = C. (B.2) 



247r2 



TT 



At 1/Q^ (the dimension two condensate is taken to be zero as implied by the sum rules in 
the S, P channels) 

(y = a{C + A^)-C Qa + j + A^A^^ = 0, (B.3) 



-16- 



= a(C + Ai)-C Qa + M^) + A^A^^ = -fl (B.4) 

At 1/Q6 

= -2a{C + A'f)iM^ + A);;) + C (^M^ (M^ + a) + ^a^j - 2CA^aW - 2A]^ Ag^ 
4 = -2aiC + A^){M^ + Ai) + C (^M^ (M^ + a) + ^a^^ - 2C7<AW - 2A^Ag) 

= i^((G;i.)'). (B.6) 



At in n^(Q2) _n^(Q2^ 

- cf = 3a(C + ^^)(m2 + Alf - 3a(C + A^f){M^ + + 6C(^:;: - Ai)A^i^ 

+ 3(A^ - Aj.)AP = -Unasimf. (B.7) 

In Eqs. (B.3)-(B.7) the following notations were introduced (the symbol i denotes m or F 
or G) 

.(3) .(3). . ar-a(a + 2M2) + ae-^n3a + 2M2) + (M2 + a)2(e^i-l)21 

Af) ^ Af )(M) = i \eB^-l)3 ^^-^0) 

In the scalar case we can write down analogous to (B.l) expansion: 

icPu-^-r (BID 

Comparing (B.ll) with OPE (2.5) and (2.6) one obtains the asymptotic sum rules. We 
write them for the 7r-out case. The 7r-in case can be easily obtained too. 



At 1/Q 

s,p 

~ 327r2 37r 

At 



■ '-^ - ' 1 + = —. (B.12) 



cf = a{C + A'a)iM' + (m^ (M^ + a) + ^a^) + AgAg)(M) 



^{{GU")^ (B.13) 
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= ^((G;1.)'). (B.14) 

At in the Us{Q^) - Up{Q^) 

cf - cf = -3a(C + A^a){M' + A^f + 3a{C + A^){M^ + - 3C{Ai - A^)Ag^ (M) 

- 3{A^a - ^G)Ag^(M) = -187ras{qqf. (B.15) 

C. Reference formulae 

The quantities A^^-^^ {j = 1,2,3, see (B.8)-(B.10)) appear due to the sums {B > 0) 



e« - r 

n=l 



E 

oo 



oo ^ 



n=l 



(e-B - 1)2 ' 



^ (e-B-l)3 • 

n=l ^ ' 

The summation of the Hnear in n part is carried out by virtue of the following asymp- 
totic representations (to be precise,once we separated out the first state, one has to make 
the shift —>■ + a in the expressions below) 

^ 1 _ ^ , f Q^ + M^ \ _ 1 f 1 L/^A#2_M 

;^^(g2 + M2 + an)2 a ) ~a\Q^~Q^\^ '2"")^ 

(m^ - aU^ + ia^) - ^ (m^ - la) [M^ - a)} ^ O (^) , (C.l) 
2^ rn2^A^2^.^,3 =^^ 2, +-1^1, , (C.2) 



„ (g2 + M2 + an)3 2a3 

n=0 ^ ' 



where 



g2 / Q2 + ^2N w 1 1 / 2 1 



3 



2g6 



Summing Eq. (C.l) and Eq. (C.3) in Eq. (C.2), one can see the absence of term at 1/g^ 
in the S, P-corrclators. The corrections contribute only to the following terms, beginning 
with l/Q^- Thus, in the scalar case the asymptotics corresponding to a dimension two 
condensate vanishes identically and we do not have a relevant sum rule. 
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D. D-wave vector mesons 



The string model does not take into account the spin of the quarks. This introduces two 
type of states depending on the relative angular momentum. In short, there exist S-wave 
and D-wave /?- mesons and this results in doubling of the p-meson trajectory [1,2]. Although 
conceptually simple, this complicates the analysis considerably. 

The two kinds of /9-mesons should enter the sum rules independently. This observation, 
generally speaking, leads to the constrain in Eq. (3.9) 

^SM + £>1 = £>), n^oc. (D.l) 
Oy ttD 

where the label D stands for the D-wave /9-mesons. It is seen from numerical estimations 
(especially for the correct ansatz ay = cld = ct) that FD{n) must be very small compared 
with Fv{n). This follows also from consideration of the matrix element {0\q^i_iq\p^). It is 
clear qualitatively that the D-wave final state is strongly suppressed compared with the S- 
one, although a concrete numerical answer will depend on how one models the wave function 
of /9°-meson. In the non-relativistic case this statement is trivial: S-wave component of 
wave function behaves as a constant and D-one does as at small distances r. Hence 
the latter tends to zero since annihilation is a point-like process. Thus, D-wave p-mesons, 
at least asymptotically, should drop out from the sum rules in contrast to what has been 
stated in [19], where this doubling was examined by means of quasi-classical string analysis. 
Let us consider the linear ansatz for the mass spectrum of D-wave vector mesons 

mjjin) = Mj) + an. (D.2) 

Direct insertion of this ansatz into sum rules (4.6) would automatically imply by the CSR 
that iJiD = M. This contradicts phenomenology [1,2]: S- and D-trajectories seem to have a 
constant splitting at any energy. In order to satisfy the sum rules of type (4.6), the D-wave 
residues should decrease at least exponentially. We propose the following parameterization 

F^{n) = aAoe"-^^", Bd > 0. (D.3) 

As these residues are exponentially small we do not consider the exponential corrections 
to the mass spectrum of D-wave /?-mesons (D.2) in our approximation. 

There arises the question as to how the chiral limit affects the behavior of the D-wave 
vector trajectory. One may think of two possibilities: i) the D-wave mesons approach the 
S-wave vector trajectory, implying asymptotic degeneration; ii) D-wave /9-mesons decouple. 
The latter variant seems to be more plausible since the former one signifies the doubling 
of states on the vector trajectory compared with the axial-vector one. In addition, the 
decoupling must take place in any case by virtue of the qualitative arguments mentioned 
above. 

E. Electromagnetic pion mass difference Am^r and chiral constants Lg, Lio- 

Given an ansatz for mass spectrum nijin) and decay constants -Fj, one can calculate the 
electromagnetic pion mass difference Amyr = — m^o., the chiral constant Liq [46] 
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(parameterizing the decay vr ei"y), and a — > vr matrix element of the electromagnetic 
penguin operator . An example of such a calculation for the LTA was provided in [16]. 
Besides, in the scalar sector we can calculate the chiral constant Lg [46] (parameterizing the 
ratio of current quark masses). We did not consider the operator Q^^^ since this quantity 
strongly depends on the cutoff and is of a limited phenomenological interest. The other 
quantities are determined by 

Am, = - ^ J° / dO'Cf {U^ld') - U\Q')) , (E.2) 

In our approach, one can directly sum over all resonances in this quantities due to the 
exact cancellation of divergent parts. Making use of notations (4. 10), (4. 13) and (D.2) one 
can write for (E.1)-(E.3) the convergent sums (we pick out the ground states, treat them 
exactly and retain only the leading and linear in exponential terms for the other ones) 



r " 



m 

oo 



rn\n)e-^-- [Aj - A^) - Ce"^-" {Aj - Q ) 



n=l 



Am. = I iC + A^) (M^ + At) In^^- 

oo 



2 / \ 

rriyiji) 



{C + A^)(m2 + A);:) ln ^'+/- + f; (Cm\n) In 

^ n=l ^ 

+ [Ce-^-- « - O + m2(n)e-^- (^^ - ^)^)] In^) } . (E.5) 



In the TT-in case 



(E.6) 



16(gg)^ 
while for the 7r-out case 

n=0 

The parameter /i in Eq. (E.5) is a normalization scale. The result does not depend on it 
as it is seen from the difference of sum rules (B.5) and (B.6). 

Accepted estimates for Liq and Lg from phenomenology are: 
^lolphen = (-5.5 ± 0.7) • 10^3 [35] ^gipj^^^ ^ (0.8 ± 0.3) • lO'^ [35]. 
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F. Numerical analysis 



In this Appendix we give an example of the meson mass spectra resulting from our work. 

The inputs general for all tables (if any) are: a = (1120 MeV)^, (qq) = -(240 MeV)^, 
iriiGt^uf) = (360 MeV)^, U = 103 MeV, = 2^, = 0.3. The units are: m(n), 
F{n), G{n) - MeV; - MeV^; Ap, Ag, BF,G,m — MeV^. 



CclS6 


lllJJUlo 


P ILrb dilQ L-Ulib LdilLh 


VA 


mviO) = 770 (769.3 ±0.8), 
mA{0) = 1200(1230 ±40), 


M = 920, B„i = 0.97, Bp = 0.72, 
AY, = -5002, j^A ^ 7702^ 

= 0.0012, A^ = -0.0031, 
Lio = -6.5 • 10-3, Am^ = 2.3 


SP 
(TT-in) 


ms{0) = 1000, 
mp(0) = 0, 

mp{l) = 1300 (1300 ±100), 
Bm = 0.97 


M = 840, Bg = 0.42, 
Ai = 5502, j^p ^ _84o2, 

A^ = -0.0009, A^ = 0.0004, 
Ls = 1.0 • 10-3 


SP 

(TT-OUt) 


msiO) = 1000, 
mp{0) = 1300(1300 ± 100), 
mp(l) = 1800(1801 ± 13), 
Bm = 0.97 


M = 1470, Bg = 1.27, 
A^ = -10802, A^ = -6902, 
A^ = 0.0213, A^ = 0.0067, 
Lg = 0.9 ■ 10-3 



Table 1: An example of parameters for the mass spectra of our work. The corresponding experi- 
mental values [2,24] (if any) are displayed in brackets. 
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2150 ^2270 ± 501 
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125 


128 


130 


msyn) 




1 nnn ^'osn-i-i 
luuu \^yoU3i iu J 




1 son n 71 


91 nn 






1 7fi 


-L 1 O 


1 78 

-L 1 O 


1 7Q 


Til 1-1 1 Tl 1 
1 1 h jJ \l I \ 




n 


1 300 (^ 300 + loo^i 


1 7fi0 (^ 801 + 1 31 


91 00 ^9070 + 351 


Gp{n) 






179 


179 


179 


m,s{n) 




1000(980±10) 


1730(1713±6) 


2120 


2420 


Gs{n) 




243 


199 


185 


181 


mp{n) 





1300 (1300 ± 100) 


1800 (1801 ± 13) 


2150(2070 ±35) 


2430 (2360 ± 30) 


Gi>{n) 




201 


186 


181 


180 



Table 2: Mass spectrum and residues for the parameter sets of Table 1. The known experimental 
values [2,24] are displayed in brackets. Not all the masses of scalar mesons are related to exper- 
imental ones since their correspondence is not well established yet because of strong mixing and 
unitarization effects. When comparing the predicted masses to the physical mass spectrum no at- 
tempt has been made to correct for the non-zero quark mass other than for the pion. Except for 
the pion the relative effect of quark masses is very small. 
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